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Abstract — This paper presents a method to improve the 
Voltage Stability of the radial distribution system by Optimal 
Capacitor placement using L-Index Sensitivity Matrix. The 
proposed methodology calculates a sensitivity matrix, which 
gives the capacitor location; it is effective than the conventional 
optimization or meta-heuristic optimization usage. The objective 
of the capacitor placement is to enhance the voltage stability. 
The proposed sensitivity matrix approach requires calculation 
of the matrix once and using it for multiple load conditions, 
against the iterative solution approach of the optimization 
methods. The characteristics of the radial distribution system is 
considered and have been used for modifications in the present 
approach. 

Index Terms — Voltage stability, Sensitivity matrix, 
Capacitor Placement. 


I. Introduction 

According to voltage and power levels, an electric power 
system as classified into a generating system, a transmission 
system and a distribution system. In many utilities, even 
though new generators have been added to augment the 
existing capacity for meeting the increase in demand, but the 
transmission capacity has not changed much due to 
constraints like rights of way and cost of tower etc. Hence, 
transmission systems are forced to carry increasingly more 
and more power over long distances using the FACTS 
devices. Voltage stability analysis is important to assess the 
criticality of the existing networks and voltage Instability is 
a very challenging problem being faced by power systems 
all over the world. 

The transmission system which is quite different, in both 
its characteristics and operation, from the distribution system. 
We cannot use same analysis approach that used in 
transmission system to distributions system. The analysis of a 
distribution system is an important area of activity, as 
distribution systems provide the vital link between the bulk 
power system and the consumers. A distribution circuit 
normally has a primary or main feeders which further rooted 
into lateral distributors and sub lateral distributors. Radial 
systems are popular because of their simple design and 
generally low cost. 
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It has been observed that the distribution system experiences 
voltage collapse under certain critical loading conditions, 
Brownell and Clarke [1] have reported the actual system 
phenomenon in which system voltage collapses periodically 
and reactive compensation is to be supplied to avoid repeated 
voltage collapse. In distribution networks, voltage instability 
is a local phenomenon and occurs near the buses which have 
high variation in loading and low-voltage profiles. So voltage 
stability analysis is very much important in distribution 
system as in transmission system for reliable and secured 
operation of the system. 

In Literature most of the work has been done on transmission 
system, voltage stability analysis [2-10], but very less work 
has been done on the voltage stability analysis in the radial 
distribution system. G.B.Jasmon and L.H.C.C Lee, [11-13] 
discussed about voltage stability analysis using various load 
flow analysis and network reduction technique and 
maximizing voltage stability by loss reduction. In [14] 
M.Moghawemi, M.O.Faruque proposed a technique to study 
the voltage collapse situation by reducing the network to a 
single line equivalent circuit and then represented through a 
7i-network. A stability index is developed to indicate the 
severity of the loading situation of the system through which 
voltage collapse is determined. An artificial Neural Network 
based technique [15] is used to predict the voltage stability 
and the learning is done is based on a new voltage stability 
index for assessment of radial distribution systems taking 
loading values P and Q as the inputs. Catastrophic theory [16] 
is proposed to identify how the equilibrium voltage point of 
the system is changed as the system parameters change during 
sudden variations in system operation. This paper explores 
the applicability of the catastrophe theory to identify a new 
VSI (voltage stability index) and finding the critical voltage 
and maximum loading for each bus in the system in terms of 
one of the elementary catastrophe. In this paper different 
indices which are used in voltage stability assessment of the 
radial distribution system is discussed. 

II. Voltage Stability Indices 

Voltage Instability is a very challenging problem being 
faced by power systems all over the world. In many utilities, 
even though new generators have been added to augment the 
existing capacity for meeting the increasing load, the 
transmission capacity has not increased correspondingly due 
to problems like rights of way, environmental concerns etc. 
Hence, transmission systems are forced to carry increasingly 
more and more power over long distances. 

An Index has to be calculated for assessing the voltage 
stability and by ranking the buses, the most prone bus for 
voltage instability is chosen. This bus becomes the location 
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for the Capacitor to be placed in the system and the 
sensitivity indices for the buses remain same as long as the 
system configuration remains the same. 

A. Line Flow Indices: 

The four line flow indices are calculated as follows. 
Line flow index for the sending end real power: 



Line flow index for receiving end real power: 

LFIRP A^[-P i+ ^Q^ ( 2 ) 

Line flow index for sending end reactive power: 

LFISQA^Q r +^P^ (3) 

Line flow index for receiving end reactive power: 

CFIRQ'A^j-Q, + ^4 2 J < 4 > 

where 

Vi :- voltages at the sending end node i 

V i+/ :- voltages at the receiving end node i +1 

r/:-resistance of the line joining the nodes i and i +1 

A7:-reactance of the line joining the nodes 

Pi :- real power flowing from node i 

Qi reactive power flowing from node i 

P, :- real power entering from node i +1 

Q r :- real power entering from node i +1 

These line flow indices are calculated for all the lines in 

the system and the hnes with a high value (close to 1) are 

considered as critical lines. The receiving end bus of the 

critical line/lines is identified as the weakest bus from the 

voltage stability perspective. 

Issues with Line flow indices: 

Consider the expressions for LFISP and LFIRP. If i, = 0 then 
LFISP = 0 and LFIRP = 0 indicating that the hnes are stable 
in the case of real power flow which could be misleading. 
This is because, the power transferred over a hne is given by 


TABLE! Line Flow Indices of 15- Bus System 


From 

To 

LFISP 

LFISQ 

LFIRP 

LFIRQ 

1 

2 

0.1075 

- 0.1111 

0.1048 

-0.1082 

2 

3 

0.0495 

-0.0505 

0.0483 

-0.0493 

3 

4 

0.0185 

-0.0186 

0.0180 

-0.0182 

4 

5 

0.0025 

-0.0025 

0.0017 

-0.0017 

2 

8 

0.0081 

-0.0081 

0.0054 

-0.0055 

8 

9 

0.0026 

-0.0026 

0.0018 

-0.0018 

2 

6 

0.0445 

-0.0454 

0.0299 

-0.0305 

6 

15 

0.0055 

-0.0055 

0.0037 

-0.0037 

6 

7 

0.0032 

-0.0032 

0.0021 

-0.0021 

3 

10 

0.0166 

-0.0168 

0.0112 

-0.0113 

10 

11 

0.0103 

-0.0103 

0.0069 

-0.0070 

11 

12 

0.0033 

-0.0033 

0.0022 

-0.0022 

4 

13 

0.0057 

-0.0057 

0.0039 

-0.0039 

4 

14 

0.0061 

-0.0062 

0.0041 

-0.0042 


B. L —index: 

In this method the Y-bus matrix of the system is split 
into rows and columns of generators and load buses. 



FLG = [YLL]-1 [YLG] (8) 




P = JLjUgjji^ (5) 

Xi 

where 8:- difference in the angles of the sending and 
receiving end voltages. 

Say the resistance of a line feeding a PQ bus ‘B’ of a 
system is zero. Let the maximum P toad at the bus B at a 
fixed value of Qi oad be P max . If the indices LFISP and 
LFIRP are calculated for the line ’L’ under these 
conditions, the values turn out to be zeros, since the 
resistance of the line is zero. Hence the line L will be 
classified as stable as per the real power loading based on 
this index whereas the system is in the verge of collapse. 


where the subscript, 

“G”:- refers to the generator buses in the system. 

“L”:- refers to the load buses in the system. 

An L-index value away from 1 and close to 0 
indicates a large voltage stability margin. The maximum 
of L-indices (L max ) of the buses to which it corresponds is 
the most critical bus. Also the summation of the squares of 
the L indices of individual buses (£L2) is used as a relative 
indicator of the overall voltage stability of the system at 
different operating conditions. 
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mode ‘V are defined as \n 


TABLE II. L- Indices of 15- Bus System 



C. Participation factor: 


It is developed as follows. Consider the load flow 
Jacobian. 


r^Li 

r H 

N ) 


UeJ 

[m 

L) 

UvJ 


Although both P and Q changes affect system conditions, 
it is possible to study the effect of reactive power 
injections on the voltage stability by setting AP (P 
constant) and deriving the Q — V sensitivities at different 
loads. Thus 


AQ = [L-MH~ l N]AV = J R AV 


(ll) 


J R = 4 A17 


where 


( 12 ) 


A : Left Eigen matrix of J R 
r| : Right Eigen matrix of J R 
E, - Eigenvalues J R 

The participation factors for the bus ‘k’ and the critical 



III. Sensitivity Matrix 


L- Index as introduced in section. 2. would predict voltage 
stability. 

To derive L- Index sensitivities Let (1,2,3 g} be no of 

generator nodes (g+1, g+2, g+3 n} be no of load nodes 

The L-index is, 

L. = 1- VF.. — (13) 

I w Vj | 

where the subscript 

“G” refers to the generator buses in the system 
“L” refers to the load buses in the system 
Squaring the equation then. 


Put, 

K=L) < 15 > 

K [> X'O.J/) (16> 

Differentiating the above equation, the new equation is: 



As the inclusion of an extra capacitor would change only 
the reactive power, we would make AP in to be zero. From 
this. 


AQ = [l-MH~ 1 n\ AV = J r AV AV = [j R \ X AQ (19 
) 

where (JR) is the reduced Jacobian. 

where AQ is a vector of change in reactive power 

injections 

AQ = [AQ g+1 AQ g+2 AQ g+3 ] ( 20 ) 
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where AV is a vector of change in voltages 

AV = [AV g+l AV g+2 AV g+3 ] 


AVj = Y j a ]i AQ i 

w = M/ 


(23) 

(24) 


As A Kj = 2 Lj A Lj 

( 26 ) 


IV. Algorithm for Capacitor Placement 

• Calculate L- Index of all the load buses as given in 

section 2. 

• Calculate the sensitivity matrix as given in section 3. 

• Choose a limit for the L- Index, let it be L-limit 

• Find the buses which is having L- Index more than 

L-limit, let the buses be [B2 B5 B7]. 

• Find A Lr ed , ALr ed is taken for the buses which are 

exceeding L-limit. 

• A Lr ed = [AL2 AL5 AL7], 

• Find the reduced sensitivity matrix (L rq ), reduced 

sensitivity matrix would relate A L rq and AQr ed where 
AQr ed =[Q2Q5Q7], 

• ALr ed = ALr d AQr ed . 

• Now perform optimization minimizing objective 

function satisfying the constraints 

AL rin < AL < 

• where 

. A L mn =L actual -L ]imit 

. A/ -_ =[ 11 -/ Vmui . 

• The optimal location of capacitor along with the Q value 

is found out. 



This is called L- Index sensitivity matrix. 

TABLE m. Comparison of a L obtained both in conventional 

AND L-INDEX SENSITIVITY 



V. Case Study and Results 

The method is carried out for finding the capacitor 
location, in 15 bus system with four times the normal load 
and the results are presented below. 

Buses exceeding L Umi are [B 4 B 5 B W B U B 12 B I3 B 14 \. 

AL red = [AL4 A L5 ALIO ALII AH 2 AH 3 AL14], 

TABLE IV. Reduced Sensitivity matrix of 15 bus experiment 
system 



AL ma;c = [0.53262 0.54971 0.55208 0.63254 0.66014 
0.57359 0.57676] 
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AL min = [0.032616 0.049709 0.052083 0.13254 0.16014 
0.073588 0.076765] 

Imposing the condition ALmin < AL < ALmax and 
running optimization program we get values of capacitors 
as given below: 


TABLE V. Capacitor Placement for Voltage Stability 
Improvement 


Bus no 

Value of Capacitor 

4 

2.3036e-011 

5 

1.8078e-011 

10 

1.6882e-011 

11 

1.5086e-010 

12 

0.38191 

13 

1.2864e-011 

14 

0.00016339 


[7] , W. Xu and Y. Mansour, “Voltage stability analysis using generic 

dynamic load models,” IEEE Trans, on Power Systems, vol. 9, pp. 
479-493, Feb. 1994. 

[8] . O. Crisan and M. Liu, “Voltage collapse prediction using an improved 

sensitivity approach”, Electric Power Systems Research, vol. 28, no. 
3, pp. 181-190, Jan. 1994. 

[9] . J. Lu, C. W. Liu, and J. S. Thorp, “New methods for computing a 

saddle-node bifurcation point for voltage stability analysis,” in IEEE 
Trans, on Power Systems, vol. 10, pp. 978-89, May 1995 

[10] T. Van Gutsem and R. Mailhot, “Validation of a fast voltage stability 
analysis method on the Hydro Quebec system”, IEEE PES Winter Power 
Meeting, 96 WM 280-8-P W RS 


The values of capacitor to be placed at bus 12 is 
0.38191MVAR and at bus 14 is 0.00016339MVAR. This 
would reduce the values of L- Index to maximum of 0.5. 
This would make the system stable. 

VI. Conclusion 

This paper emphasis on the analyzing and improvement of 
voltage stability of radial power systems. The various indices 
which have been proposed for assessing voltage stability are 
studied. Patterns of voltage stability indices are been 
observed and their dependence on the electrical distance has 
been studied. An algorithm is proposed for optimal placement 
of the capacitors to improve the voltage stability to the 
required level based on the new sensitive index proposed. The 
effect of reactive power injections at a bus on the entire power 
system in consideration is studied. The algorithms which have 
been proposed for radial power systems can also be suitably 
modified and can be applied to meshed systems. 
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